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ABSTRACT 

Symmetry adaptation techniques are applied to the determination 
of the intensity of two-photon (intra- and inter-configurational) 
transitions for transition ions in finite symmetry environments. 

1. INTRODUCTION 

Symmetry adaptation techniques, developed in the spirit of Refs. [1,2], for a 
chain of groups 0(3) D G are applied here to two-photon spectroscopy of transition 
ions of configuration n£ N (£ = d for transition metal ions and £ — f for lanthanide 
or actinide ions) in surroundings of finite symmetry G. More precisely, we show in 
this lecture how Wigner-Racah calculus for a chain 0(3) D G (in terms of simple 
or double groups) can be combined with models 3-20 -* based on second- plus third- 
order mechanisms in order to isolate the polarization dependence from the intensity 
of two-photon transitions for an n£ N ion in a molecular or solid-state environment 
with symmetry G. (For classification and symmetry-breaking purposes, the group 
G may be replaced by a chain of subgroups of 0(3), the relevant symmetry group 
being one of the groups of the chain.) 

The role of symmetries in two-photon spectroscopy of partly-filled shell ions 
in finite symmetry is touched upon in Refs. [4-6,12-14,16,17]. In Refs. [4,5], the 
information arising from symmetry is handled mainly in a qualitative way. More 
quantitative results can be found in Refs. [6,12-14,16,17]. In Ref. [12], the accent is 
put on the transition matrix elements between initial and final state vectors while 
emphasis is on the intensity strength in Refs. [14,16] and in the present lecture. 
This lecture constitutes a complement to the material presented in Refs. [14,16]. 

Two distinct cases are studied in this work. The case of (n£ N — > n£ N , e.g., 
3^ — > 3^ and Af N — > 4/^) intra-configurational two-photon transitions, which 
are parity allowed, is worked out in section 2 and the one of (n£ N — > n£ N ~ l n'£' 
with £ + £' odd, e.g., 3^ — > 3d N ~ 1 4p and Af N — > 4f N ~ 1 bd) inter-configurational 
two-photon transitions, which are parity forbidden, is examined in section 3. 

2. INTRA-CONFIGURATIONAL TWO-PHOTON TRANSITIONS 
2.1. Preliminaries 

We know that the electronic transition matrix element Mi_>f between an 
initial state i and a final state / is, in the framework of the electric dipolar ap- 



proximation, given by 



M ^f = {f\D-^\v) ( V |£.£|t) {v\D.£2\i) (1) 

v 1 v 2 

The two summations in (1) have to be extended over all the intermediate states 
v having a parity different from the one of the states i and /. Furthermore, we 
have A a = Tiui\ — E v , where E v is the energy of the state v with respect to that 
of the state i and fiu x the energy of the photon no. A. (For Raman scattering, the 
sign of hu2 has to be changed.) The quantity D.£ x in (1) stands for the scalar 
product of the electric dipolar moment operator D for the N electrons and the 
unit polarization vector £ x for the photon no. A. (We use single-mode excitations, 
of energy Tiu x , wave-vector k x and polarization £ x , for the radiation field.) The 
two photon beams can be polarized either circularly with 

(€ x ) q = -S(q,-l) if Sx = e +1 while (S x ) q = -S(q, +1) if £ x = (2) 
or linearly with 

(£ A ) = cos6» A (£a)±i = =F^= sin^A exp(±i(f X ) (3) 

In equations (2) and (3), we use the components {£ x ) q = £ x .e q (with q = — 1, 0, 1 
and A = 1, 2) in the standard spherical basis (e_i, eb, e+i). In the case of a linear 
polarization, the angles (9 X , <p x ) are the polar angles of the polarization vector £ x 
(A = 1,2) with respect to the crystallographic c-axis. For two-photon absorption, 
only one sum occurs in (1) when the two photons are identical. 

Equation (1) can be derived from the time-dependent perturbation theory 
and goes back to the work of Goppert-Mayer 18 ' 19 ) (see the lecture by J.C. Gacon 
in these proceedings). It is also possible to derive it, in an elegant way, from the 
method of the resolvent operator 20 - 1 . 
2.2. State vectors 

The initial state i with symmetry V is characterized by the state vectors ^7) 
where 7 (7 = 1, 2, • • •, dim T) is a multiplicity label to be used if the dimension 
dim T of the irreducible representation class (IRC) V of the group G is greater 
than 1. The state vector ^7) is taken in the form 

|zr 7 ) = \rd N iT-f) = \n£ N aSLJaT-f) c(aSLJaT;i) (4) 

aSLJa 

in terms of the 0(3) D G symmetry adapted state vectors 1 ' 2 ) 

J 

\nl N aSLJaT-f) = \n£ N aSLJM) (JM|Jar 7 ) (5) 

M=-J 



The coefficients ( JM\ JaTj) in (5) are reduction coefficients to pass from the chain 
0(3) D 0(2) characterizing the {JM} scheme to the chain 0(3) D G character- 
izing the {JaF'j} scheme ; they depend on the group G with a certain degree 
of freedom emphasized by the branching multiplicity label a to be used when V 
occurs several times in the IRC (J) of 0(3). In contradistinction, the coefficients 
c(aSLJaT; i) in (4) depend on the Hamiltonian employed for obtaining the initial 
state i. Similarly, for the final state / with symmetry r', we have the state vectors 

\fY'i) = \n£ N fY'i) = Yl \m N a'S'L'J'a'Y'i)c{a'S'L'J'a'Y'-f) (6) 

a'S'L'J'a' 

in terms of 0(3) D G symmetry adapted state vectors. The only good quantum 
numbers for the initial and final state vectors are 1?7 and T'j', respectively. Al- 
though, the state vectors \iT-f) and |/r" 7 ') are developed in a weak-field basis, it 
is to be noted that the intensity calculation to be conducted in what follows is 
valid for any strength (weak, intermediate or strong) of the crystalline field. 
2.3. Transition matrix element 

By using a quasi-closure approximation, it can be shown that the transition 
matrix element M^f between the state vectors |£T 7 ) and \fT /r j') is given by 

M^ f = M i(r7W(r , 70 = (/rVl^e//|<r 7 ) (7) 

where H e ff is an effective operator 7 ' 8 ). This operator may be written as 12 ) 

H 'ff= E E C[{k s k L )k] fas^K W^) fc ) (8) 

fc=0,l,2 k s k L 

In equation (8), W^ sfci ^ fc is an electronic double tensor of spin rank ks, orbital 
rank k^ and total rank k. The information on the polarization of the two photons 
is contained in the tensor product {£1 £2} of rank k = 0, 1 or 2. The right- 
hand side of (8) is a development in terms of scalar products ( . ) with expansion 
coefficients C[(kski)k]. These coefficients depend on the ground configuration 
nl N and on the configurations n£ N ~ 1 n '£' and/or n'£' u +1 n£ N+1 , with £ + £' odd, 
from which the states v arise. 

Only the contributions (ks = 0, ki, = 1, k = 1) and (ks = = 2, k = 

2) correspond to the standard theory originally developed by Axe 3 ). The other 
contributions (ks 7^ 0, kL,k), which may include (ks = 1, ki, = 1, k = 0) and 
(ks = = l,k = 2), correspond either to mechanisms introduced by various 
authors 7-11 ) or to phenomenological contributions 12 -*. The contributions (ks = 
0, kL = k, k) and (ks 7^ 0, k) are often referred to as second-order and third- 
order mechanisms, respectively. It is in principle possible to find an expression for 
the parameters C [(kski,) k]. Among the various contributions (ks 7^ 0, k]^, k), the 
contribution (ks = 1, &l = 1, k = 0) arises from the spin-orbit interaction within 
the configuration n£ N ~ l n'£' as was shown for lanthanide ions 7 ' 8 ). 



The transition matrix element (7) is easily calculated by means of Wigner- 
Racah calculus for the chain 0(3) D G. As a result, we have 12 ) 

M i( r 7 )->/(r'7') = £ £ c(a'S'L'j'a'T'-J)* c(aSLJaT;i) 

a'S'L'J'a' aSLJa 

(_)ks+k L -k C[(k s k L )k] (n£ N aSLJ\\W( ksk ^ k \\nl N a'S'L'J'y 
E/( a r 7 a'Y'i a"Y"i' ) ^ £ ^a"T"i" 



ksk L k 



a"T"j" 

where the / symbol denotes an 0(3) D G symmetry adapted coupling coefficient 
defined by 1 ' 2 -* 

(J J' k \ v (U-m( J k J' 

J I ar 7 aTV a'T'V' ) 2^ ^ > \ -M q M' 

(JM|JaI»* (J'M'I J'aTV) (fc#a"r" 7 ") 

Equation (9) follows by developing (7) with the help of (4), (6) and (8). 
2.4. Intensity formula 

The quantity of interest for a comparison between theory and experiment is 
the intensity S , ^r)^/(r / ) of the two-photon transition between the initial state i 
and the final state /. This intensity is given by 

S T ^v = S^rWT') = £ l M *(r7W(r'7')| 2 ( 10 ) 

77' 

By introducing (9) into (10) and using the factorization property 1 -* for the / co- 
efficients as well as the orthonormality-completeness property 1 ) for the Clebsch- 
Gordan coefficients of the group O, we obtain the compact expression 

= EEE I[kirsT";TT'} £ {£i*2}$ v , (Vi£ 2 }S v ,)* (11) 

k,£ r,s T" 7 " 

In equation (11), the parameter / reads 



(12) 



I[klrsT"-TT>] = [V"]- 1 £ £ £ £ Yk(J'a'V, JaT) 

J'a' Ja J'a> Ja 

Y £ (J'a'T', JaT)* £ ( J ' a ' r ' + krT"\Ja(3T) {J'a'V + lsT"\Ja[3T)* 

H 

where Yk is defined by 

Y k (J'a'T',JaT) = [J}~ 1 / 2 £ £ £ c(a'S'L' J'a'V; /)* c(aSLJaT; i) 

a'S'L' aSLk s k L (13) 

C[{k s k L )k] (-) k s+k L -k (ne N aSLJ\\W( kskL)k \\ne N a'S'L , J'y 



and Yi by a relation similar to (13). In (12) the ( + | ) coefficients stand for 
isoscalar factors of the chain 0(3) D G and the labels (3 are internal multiplicity 
labels to be used for those Kronecker products which are not multiplicity-free 1 -'. 
2.5. Properties and rules 

The / parameters in (11) can be calculated in an ab initio way or can be 
considered as phenomenological parameters. In both approaches, the following 
properties and rules are of central importance. 

Property 1. In the general case, we have the (hermitean) property 

I[eksrT";TT']* = I[k£rsY"- FT'] (14) 

which ensures that Sr->r' * s a rea l quantity. 

Property 2. In the case where the group G is multiplicity- free, we have the 
factorization formula 

I[k£rsT";TT'] = X [krT"; TV'} x[tsT"; TV']* (15) 

where the function \ is defined through 

x[£xr" ; rr'] = [r"]- 1/2 [r] 1 / 2 ^ ^ Y k (j'a'r', j a r) ( j'a'V + krY"\jaY) 

J'a' Ja 

(In a less restrictive sense, equation (15) is valid when the Kronecker product 
r'* <g> T, of the complex conjugate IRC of V by the IRC T, is multiplicity- free.) 

The number of independent parameters / in the expansion (11) can be a 
priori determined from the two following selection rules used in conjunction with 
Properties 1 and 2. 

Rule 1. In order to have S'rwr' 7^ 0, it is necessary that 

r" cr'*®r r" c (k g ) r" c (e g ) (16) 

where (k g ) and (£ g ) are gerade IRC's of the group 0(3) associated to the integers 
k and £, respectively. 

Rule 2. The sum over k and £ in the intensity formula (11) is partially 
controlled by the selection rule 

E\ 7^ £2 '■ k, £ = 1, 2 for 2nd-order or fc, £ = 0, 1, 2 for 2nd-order + 3rd-order 

or 

E\ = £2 : k, £ = 2 for 2nd-order or k, £ = 0, 2 for 2nd-order + 3rd-order 

according to as the two photons have different or the same polarization. (Note 
that the situation E\ = £2 surely occurs for identical photons but may also occur 
for non-identical photons.) 



2.6. Discussion 

For most of the cases of interest, there is no summation on r and s, two 
branching multiplicity labels of type a, in the intensity formula (11). (In other 
words, the frequency of r" in (k g ) and (£ g ) is rarely greater than 1.) The group- 
theoretical selection rules (16) impose strong limitations on the summation over 
r" in (11) once V and V are fixed and the range of values of k and £ is chosen. 

The number of independent intensity parameters / in the formula (11) de- 
pends on : (i) the nature of the photons, cf. Rule 2 ; (ii) the group G, cf. Rule 
1 ; (iii) the conjugation property (14), cf. Property 1 ; (iv) the use of ks = 
(second-order mechanisms) or ks = and ks ^ (second- plus third-order mech- 
anisms), cf. Rule 2 ; (v) the (weak-, intermediate- or strong-field) state vectors 
used in conjunction with equations (12) and (13). 

Points (i)-(iii) depend on external physical conditions. On the other hand, 
points (iv) and (v) are model-dependent. In particular, in the case where the J- 
mixing, cf. point (v), can be neglected, a situation often of interest for lanthanide 
ions, the summations on k and £ in (11) are further reduced by the triangular rule 
| J — J' | < k, I < J + J', where J and J' are the total angular quantum numbers 
for the initial and final states, respectively. Similar restrictions apply to ks and 
kh in (13) if the S- and L-mixing are neglected. 

The computation, via equations (12) and (13), of the I parameters generally 
is a difficult task. Therefore, they may be considered, at least in a first step, 
as phenomenological parameters. In this respect, equations (12) and (13) should 
serve as a guide for reducing the number of I parameters. 

Once the number of independent parameters / in the formula (11) has been 
determined, we can obtain the polarization dependence of the intensity strength 
Sr^v by calculating the tensor products {Si £2}a»r» 7 » (with K = k,£ and a" = 
r, s) occurring in (11). For this purpose, we use the development 

{£i£ 2 }S„ 7 „ = J2 {Wg {KQ\Ka"Y"i') (17) 

Q=-K 

in terms of the spherical components {Si £2)0^, the coefficients in (17) being re- 
duction coefficients for the chain 0(3) D G. Then, we use in turn the development 

{^Mif = R K -W 2 t t (L - Q 

m= — 1 m' = — 1 ' 

in terms of the spherical components (£\) q defined by (2) or (3) for circular or 
linear polarization, respectively. 

2.7. Illustration 

As a pedagological example, let us consider the case of the two-photon ab- 
sorption (intra-configurational) transition 7 Fq — > 5 _Dq for the configuration 4/ 6 in 
tetragonal symmetry (with G = C^ v or D 2 d)- In this case, we have r = A\ for the 



initial state and V = At for the final state, whatever the strength of the crystalline 
field is. Rule 1 then yields V" = At and, consequently, there is no sum on the label 
7" in (11). Furthermore, there is no sum on the branching multiplicity labels r and 
s in the intensity formula (11). Let us consider an experimental situation where 
the two photons are identical (one-color beam arrangement) so that £\ = £ 2 = £• 
We continue with a model characterized (in addition to the symmetry C^ v or Did) 
by the use of second- plus third-order mechanisms. Thus, according to Rule 2, the 
indices k, t in (11) may assume the values and 2. By introducing the abbreviation 
I{kl) = I[k£rsF" = At; T = A\ V = Ai], we are left with 3 a priori independent 
parameters (cf. Property 1), viz., 7(00), 7(02) = 7(20)* and 7(22). Taking the 
wave-vector k of the two photons parallel to the crystallographic c-axis, we have 

— 1 3 ens 2 ft — 1 
{££}% ee {££}* = or and {££} 2 Al ee {££}l = or (18) 

for linear or circular polarization, respectively. Finally, the intensity strength 
Sai^Ai is easily obtained by introducing (18) into (11). This leads to 

S Al ^A 1 = r 2 + s(3cos 2 #-l) + t 2 (3cos 2 fl- l) 2 or (19) 

according to whether as the polarization is linear or circular, respectively. Here, 
the real parameter s and the two non-negative parameters r 2 and t 2 are defined 
by 

2 _ £M - -A02)-/(02)* ,2 _ £(22) 

3 3^2 6 

In the intensity formula (19), the two first terms (in r 2 and s) arise from third- 
order mechanisms while the third one (in t 2 ) comes from second-order mechanisms. 
Only the scalar term (in r 2 ) contributes to S Al ^A 1 in the absence of J-mixing. 

Equation (19) has been applied 13 ) to the cases of Sm 2 +:BaClF, Sm 2 +:SrClF 
and Eu 3+ :LuP04. Neither the conventional second-order term (in t 2 ) nor the 
scalar third-order term (in r 2 ) are sufficient to reproduce the experimental data. 
Indeed, on the basis of fitting procedures 13 ^ , a good agreement between theory and 
experiment requires in these cases that the three terms (in r 2 , t, and s 2 ) contribute 
to the intensity strength (19). As a conclusion, the model inherent to equation 
(19), the ingredients of this model being : (symmetry, second- plus third-order 
mechanisms, J-mixing), is appropriate to the 7 Fq — > 5 Do two-photon transition 
for Sm 2 +:BaClF, Sm 2 +:SrClF and Eu 3+ :LuP0 4 . 

Similar models have been applied to the two-photon transitions 13 ^ 7 7o(r = 
At) -> 5 D 2 (T' = At,Bt,B 2 ,E) for the tetragonal compound Sm 2+ :BaClF and 
to the two-photon transitions 17 ) 3 A 2 (T = T 2 ) -> 3 T 2 (r' = A 2 ,E,T 1 ,T 2 ) for the 
cubical compound Ni 2+ :MgO. In both cases, we have found that second-order 
mechanisms are sufficient to describe the transitions. More precisely, the model : 
(symmetry, second-order mechanisms, S- and L-mixing but no J-mixing) works 



for Sm 2+ :BaClF while the model : (symmetry, second-order mechanisms, S-, L- 
and J-mixing) works for Ni 2+ :MgO. 

3. INTER-CONFIGURATIONAL TWO-PHOTON TRANSITIONS 
3.1. Sketch of the theory 

We now consider two-photon transitions between Stark levels arising from the 
configurations n£ N and n£ N ~ l n'£' of opposite parities (£ + £' odd). For the sake of 
simplicity, we deal here with identical photons. The initial (i.e., ^7)) and final 
(i.e., l/r'7')) state vectors are taken in the form 

\n£ N iTj)= \n£ N aSLJaT-f) c(aSLJaT;i) 

otSLJa (or\\ 

l^-^r/r'V) = Yl W^rit'ctS'L'J'a'T'i) c{a'S'L'J'a'T'- f) 

a'S'L'J'a' 

to be compared with equations (4) and (6). 

It is clear that the transition matrix element 

M i(r7W(r , 7 ,) = X [f v 'i\D.E\vY vlv ^ (vT vlv \D.£\irJ) 

V 

is identically zero. In order to obtain Mj(r 7 )^/(r'7') 7^ 0, it is necessary to pollute 
(20), as well as the intermediate state vectors, with state vectors of the type 
\n£ N ~ 1 n'£'xT'Y) and \n£ N xTj), respectively. This may be achieved by using 
first-order time-independent perturbation theory where the polluting agent is the 
crystal-field potential H3 of odd order, which is static or dynamic according to as 
the group G does not or does have a center of inversion. Hence, we produce state 
vectors noted \n£ N iTj > and \n£ N ~ 1 n' '£' ' fT '7' > from which we can calculate, 
in a 2nd-order time-dependent plus lst-order time-independent scheme, a non- 
vanishing transition matrix element 

M i( r 7H/( r' 7 ') = X < f T 'y\D-£\vT vlv >< vT^lD.SliT-f > 

V 

Then, we apply a quasi-closure approximation both for the initial, intermediate, 
and final state vectors and the transition matrix element. This approximation can 
be summarized by E(n'£') —E(n£) = 2hu. We thus obtain a closed form formula 
for Mj(r 7 )^j( r / 7 /) (see Ref. [20]). 

At this stage, it should be mentioned that the so-obtained formula is equiva- 
lent to that we would obtain, within the just mentioned approximation, by using 
third-order mechanisms described by 

M i(r7W(r , 7 ,) = XT-T xt-t (/rYI^-^kiri7i)(^ri7il^-^b2r 2 72)x 

x (u 2 r 2 72|#3|ir7) + term [D. £ \H 3 \D. £] + term [H 3 \D.£\D.S] 



where the initial, intermediate and final state vectors are non-polluted. 

By following the same line of reasoning as in the case of intra-configurational 
transitions, we are left with the intensity formula 20 ^ 

S r ^=Re[ EE^^'^E^^Y 

k,£=0,2 r,s T" ■," ^ 

+ E EE^[^r'' ; rr^{55}« 7 „ {ee}^„] 

fc=0,2 r,s V" 7" 

which parallels the formula (11). A detailed expression of the intensity parameters 
I\ and I2 will be found in the thesis by Daoud 20 ) and in forthcoming papers. 
3.2. Illustration 

Let us consider the case of the configuration 4f in tetragonal symmetry with 
G = C^ v and examine the two-photon transitions between the Stark levels of the 
shells 4f and 5c? (i.e., n£ = 4/, N = 1, n't' = hd). There are four possible 
transitions since the initial and final states may have the symmetries T 6 and T 7 . 
For a linear polarization, the application of the intensity formula (21) leads to 

Sr 6 — r 7 = f vr 2 + g n s + h 7r 4 + % 7r 5 , Sr 7 — r> = a' + b' 7Ti + c irj + d' 7r 2 + e tt 3 
Sr 7 ^r 6 = + g' 7T3 + ^'tt4 + i'n 5 , S Tfi ^r e = a + 6tti + ctt 2 + dir 2 + eTv 3 
where the angular functions %i (i = 1,2, 3, 4, 5) are defined by 
7Ti = 3 cos 2 — 1, 7r 2 = sin 2 20, 7r 3 = 7T2Cos2(/?, 7r 4 = sin^ 9 cos 2 2y?, 7T5 = sin 4 — ^4 

The various parameters a, • • • , i and a', • • • , i' are simple functions 20 - 1 of the inten- 
sity parameters I\ and I2 occurring in (21). 
4. CLOSING REMARKS 

We have shown how 0(3) D G symmetry adaptation allows to derive intensity 
formulas for intra- and inter-configurational two-photon transitions for ions in 
molecular or solid-state environments. In particular, the number of independent 
parameters required for describing the polarization dependence of the transitions 
is determined by an ensemble of properties and rules which combine symmetry 
and physical considerations. The main results of this paper are formulas (11) and 
(21) for intra- and inter-configurational transitions, respectively. 

The polarization factors in (11) and (21) are under the control of the ex- 
perimentalist. Both formulas depend on expansion coefficients c(aSLJaT;i) and 
c(ot'S'L'J'aT'; /). These coefficients (model dependent) can be obtained by opti- 
mizing Hamiltonians, for the ion in its environment, involving at least Coulomb, 
spin-orbit and crystal-field interactions ; the introduction of more sophisticated in- 
teractions may be useful to take covalency effects into account 2 ^ . Alternatively, the 
expansion coefficients can be considered as free parameters. Furthermore, in (11) 
and (21) we have reduced matrix elements (configuration dependent), isoscalar fac- 
tors for 0(3) D G (group theory dependent), and C[(kski,)k] parameters (mecha- 
nisms dependent). As a conclusion, there are three ways to deal with the intensity 



parameters in (11) and (21) : they can be calculated from first principles, or 
considered as phenomenological parameters or determined in a mixed approach. 

The author is very grateful to the Organizing Committee of the School for 
inviting him to give this lecture. 
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